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BIRTHDAY INEQUALITIES, REPULSION, AND HARD SPHERES 


WILL PERKINS 


Abstract. We study a birthday inequality in random geometric graphs: the probability 
of the empty graph is upper bounded by the product of the probabilities that each edge 
is absent. We show the birthday inequality holds at low densities, but does not hold in 
general. We give three different applications of the birthday inequality in statistical physics 
and combinatorics: we prove lower bounds on the free energy of the hard sphere model and 
upper bounds on the number of independent sets and matchings of a given size in d-regular 
graphs. 

The birthday inequality is implied by a repulsion inequality, the expected volume of the 
union of spheres of radius r around n randomly placed centers increases if we condition on 
the event that the centers are at pairwise distance greater than r. Surprisingly we show that 
the repulsion inequality is not true in general, and in particular that it fails in 24-dimensional 
Euclidean space: conditioning on the pairwise repulsion of centers of 24-dimensional spheres 
can decrease the expected volume of their union. 


1. Introduction 

How many people must be in a room so that the chance at least two share a birthday is 
at least 1/2? This is the ‘Birthday Problem’, and the answer is that 23 people is enough 
(assuming that the birthdays are independently and identically distributed). 

Our starting point is an elementary inequality, which we will call the birthday inequality. 

Proposition 1 (The Birthday Inequality). Suppose n people have birthdays each chosen 
independently and uniformly at random from m possible birthdays. Let E n be the event that 
no two people share a birthday, and p = 1/m the probability that two given people share a 
birthday. Then 

Pr [E n \ < (1 -p)(a) . 

Proof. Let Ey be the event that there are no shared birthdays among the first k people. Let 
14 be the fraction of birthdays covered by the first k people. Then 

Pr[E n ] = E[1 - V n -i\E n _i] • Pr[E n _i] . 

We assume inductively that Pr[£/] < (1 — p)^ 2 \ and note that 

E[l-V k \E k \ = l--<(l--) =(1 -p) k 
m \ m J 
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which shows that Pr[i^ +1 ] < (1 — p)( 2 ) for all k > 1. □ 

We are interested in geometric birthday inequalities, and in particular in settings relevant 
to two models from statistical physics: the hard sphere model and the hard-core lattice gas 
model. In these models, particles are placed at random in a metric space X equipped with 
a probability measure p (e.g. the unit cube or a subset of the d-dimensional integer lattice 
with uniform measure) conditioned on all pairwise distances between particles being larger 
than some threshold r. In this setting the birthday inequality supposes an upper bound 
on the probability that no two particles are within distance r when n particles are placed 
independently at random according to p. 

Definition 1. Let Ah, X 2 ,... X n be independently sampled points from a space X according 
to the distribution p. Then a birthday inequality holds if 

(1) Pr[E n ] < (1 - p)® 

where E n is the event {Ai <i < j< n d(Xi, Xj) > r} and p := Pr[d(Ai, X 2 ) < r]. 

The quantity on the right is what the probability on the left would be if all pairwise inter¬ 
actions were independent, and so the birthday inequality is a statement about correlations 
of these events. 

A related inequality is the following repidsion inequality. 

Definition 2. In the setting above, the repidsion inequality holds if 

(2) E[V k \E k ]>K[V k ] 

where is the event that the centers X {, X 2 , . . . X^ are at pairwise distance greater than r, 
14 is the volume fraction of X covered by the union of the closed balls of radius r around 
X \,... Xk , and the expectations are taken over choosing X \,... X^ independently at random 
in X according to p. 

The repulsion inequality states that conditioning on the event that the centers of randomly 
placed balls of radius r are at pairwise distance greater than r does not decrease the expected 
volume of their union (as compared to the unconditional expectation). The repulsion inequal¬ 
ity has the flavor of a probabilistic version of the Kneser-Poulsen conjecture ushu]: moving 
a set of spheres in Euclidean space so that all pairwise distances between their centers do 
not decrease cannot decrease the volume of the union of the spheres. This was proved in two 
dimensions by Bezdek and Connelly pQ, but is open in higher dimensions. 

While the repulsion inequality seems intuitively obvious, we show in Corollary [2] that it is 
not always true; in particular it fails in dimension 24. 

As in the proof of Proposition [1] the birthday inequality on n points is implied if the 
repulsion inequality holds for all k between 1 and n—1. We write Pi'[E n ] = (1—E[V4-i|^'n-i])‘ 
Pr[E n _i] and continue inductively. The unconditional expectation satisfies E[V&] = 1 — (1 — 
p) k , and so if the repulsion inequality holds for all 1 < k < n — 1, the birthday inequality 
holds. 

We will show that at sufficiently low particle densities, the repulsion inequality holds in 
both the hard sphere and hard-core models. This leads to bounds on the free energy in 
both models via the birthday inequality. However, we will also show that at sufficiently high 
densities the birthday inequality can fail. We conclude by conjecturing that the failure of the 
repulsion inequality can be used to indicate the fluid/solid phase transition. 
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2. Hard spheres 

The hard sphere model is a model of particles as randomly positioned non-overlapping 
spheres: a random sphere packing. There are no forces in the model besides the hard 
constraint that two spheres cannot overlap. We define the hard sphere model on T d , the 
d-dimensional unit torus. 

Definition 3. The hard sphere model H d {n, r) consists of a uniformly random configuration 
of n spheres of radius r/2 in T d , conditioned on the event that the centers of the n spheres 
are at pairwise distance greater than r. 

An important quantity in statistical physics is the partition function: 

Definition 4. The partition function, Z d (n,r), of the hard sphere model on T d is defined as: 

(3) Z d (n, r)= ■■■ l En dxi--- dx n 

JT d Jr d 

where E n is the event that d(xi,xj) > r for all 1 < i < j < n. 

We define the density a of H d (n, r) as the fraction of volume of T d occupied by the spheres 
of radius r/2 around the n centers, i.e., a = n(r/2) d v d where v d is the volume of the unit 
ball in M. d . As a is the density of the random sphere packing given by H d (n,r), it must lie 
between 0 and the maximum sphere packing density in d dimensions. 

Definition 5. The free energy of the hard sphere model at density a is: 

(4) F d {a) = - lim - log Z d (n,r n (a)) 

n—>oo Ti 

where r n (a) = 2 (a/(nv d )) 1 / d . 

Physicists believe that the hard sphere model in dimension d > 2 undergoes a fluid / 
solid phase transition as the density of the spheres increases: at low densities configurations 
show no long-range order, while after the phase transition long-range order emerges. For 
an introduction to the hard sphere model see m and the references therein. In dimension 
d = 1 there is no phase transition and the model is solved; that is, an explicit expression 
for the free energy is known, and it has no non-analytic points [21] . Mathematicians have 
proved rigorous lower bounds on the density at which a Markov chain to sample from the 
model mixes rapidly DSUim See [ISU2] for a discussion of mathematical proofs of phase 
transitions in continuous hard-core models, and in the second a proof of a phase transition 
in a system with zipper-like molecules. 

We define the model with spheres of radius r/2 because it will be convenient to view the 
hard sphere model from the perspective of the random geometric graph G d (n,r ): n points 
placed uniformly and independently at random in T d with an edge placed between pairs of 
points at distance at most r. 

The following proposition relates the hard sphere model to the random geometric graph, 
and follows immediately from Definition [3l 

Proposition 2. 

Z d (n,r) = Pr [G d (n,r) is empty}. 
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In what follows we parameterize both the hard sphere model and the random geometric 
graph by p := Vdr d , the probability that two uniformly random points in T d are at distance 
at most r. Abusing notation we will write Gd(n,p) for Gd(n,r(p )) where Vdr{p) d = p. This 
parameterization gives some intuition for the birthday inequality: if G(n,p) is the Erdos- 
Renyi random graph on n vertices (every edge present independently with probability p), 
then the birthday inequality is Pr [Gd(n,p) is empty] < Pr [G(n,p) is empty]. In fact, if we 
fix n and p and let d —> oo, then Pr [Gd(n,p) is empty] —>• Pr \G{n,p) is empty] (Theorem 2 in 
[6] for the random geometric graph defined on the surface of the d-dimensional unit sphere). 

The birthday inequality holds in dimension 1 for all values of p: 

Proposition 3. For all p € [0,1], 

Pr [Gi(n,p) is empty} < (1 — p)^ 2 ^ . 

Proof. For p > 2/n, we are beyond the maximum packing density on the circle, and so 
Pr[Gi(n,p) is empty] = 0 and the inequality holds. For p < 2/n, we can write the left-hand 
side explicitly: Pr [G\{n,p) is empty] = (1 — np/ 2) n_1 ; then it is a calculus exercise to show 
that (1 — np/2) n ~ l < (1 — p)(z) for p < 2/n. □ 

Our first main result of this section is to show that in any dimension, at a low enough 
density the birthday inequality holds. We do this via the repulsion inequality (f2|) . 

Theorem 1. For the the hard sphere model on T d , for densities a < 2~ 2 ~ 5d the repulsion 
inequality holds. 

Theorem[l]and the birthday inequality immediately imply a lower bound on the free energy 
of the hard sphere model at sufficiently low densities, which to the best of our knowledge is 
new. 

Corollary 1. For a <2 2 3d , Fd(a) > 2 d l a. 

Proof of Theorem Q We first define some notation. For a collection of k centers in T d , let 
14 be the volume of points in T d at distance at most r from one of the k centers, i.e. the 
volume of the union of balls of radius r around the centers. Let E & be the event that the k 
centers are at pairwise distance greater than r. As always, we have p = Vdr d , and we assume 
a < 2~ 2 ~ 3d , i.e. p < 4~ 1_rf /n. Our goal is to prove the repulsion inequality E[T4|F4] > E[T4] 
where the randomness is in placing each center uniformly and independently at random in 
T d . 

We will prove the following estimate for all 1 < k < n — 1: 

( 5 ) nVtm > kp .Qplz£F. 

To complete the proof of the Theorem from ([5]) we use inclusion/exclusion to bound E[I4] = 
1 — (1 — p) k < kp — (g )p 2 + {^)p 3 , and so 

which is non-negative when p < 4 ~ 1 ~ d /n. 
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To prove © we use inclusion/exclusion and linearity of expectation to get the lower bound 
m k \E k ] > kp-^K[V(i,j)\E k \ = kp- QW V(l,2)\E k \ , 

i<j ' ' 


where V(i,j ) is the volume covered by the balls of radius r around both centers i and j, i.e. 
their overlap volume. 

Now let x be a fixed point in T d (say the origin), A^’ 2 the event that x is covered by the 
balls of radius r around both centers 1 and 2, and E 2 the event that centers 1 and 2 are at 
distance greater r. Then we have 


E[V(l,2)\E k ] 


Pr[4’ 2 | E k ] 

Pr[4’ 2 n E k ] 

Pr [E k ] 

Pr[4’ 2 n e 2 ] ■ Pr[S fc |4’ 2 n e 2 ] 


Pr [E 2 ] • Pv[E k \E 2 ] 


Pr[4’ 2 !^] 


Pr[E k \Al’ 2 n E 2 ] 
?r[E k \E 2 ] 


First note that Px[E k \A^ 2 n E 2 \ < Pr{E k _ 2 ]. Next, write 


Pr[E k \E 2 ] 


_ Pr [E k _ 2 ] Pi[E k \E k _ 2 } 

Pr [E 2 ] 

> Pr[£ t - 2 ] (1 — ( *~ 2 M±EtzM 

1 — p 


> Pr [E k _ 2 \ 


(1 — kp ) 2 
1 ~P 


where we have used the inequalities Pv[E k ^i\E k _ 2 \ > 1 — (k — 2 )p and Pr[E k \E k _i] > 1 — {k — 
1 )p which follow from the union bound: the volume of the union of balls of radius r around 
k — 2 centers is at most (k — 2 )p. 

This gives 

Pr[E k \Al’ 2 n e 2 \ 1-p 

Pr[E k \E 2 ] - (1 -kpf 


Finally we upper bound Pr[^’ 2 |Fl 2 ]- We write 


p 2 = Pr[4’ 2 ] =pPr[A 1 /\E 2 ] + (1 - p) Pr[4’ 2 |£ 2 ]. 


The probability that three given points form a triangle in the random geometric graph with 
connection radius r is p ■ Prf^’ 2 )^]- A lower bound for the probability of forming a triangle 
is the probability that the first two points fall in a ball of radius r j 2 around the third, which 
has probability p 2 A^ d . Putting this together we have 


Pt[A]: 2 \E 2 } < 


p 2 — p 2 4 d 

1 — p 
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and 

i^ 1 ’ 2 )' £ ‘]^ 2 (WW' 

which gives Q. 

□ 

Our next result is that the birthday inequality does not hold in general. We show this in 
dimension 24, using the fact that there is a sphere packing of particularly high density. 

Theorem 2. In dimension 24, the birthday inequality fails for large enough n at densities 
a € ((.79) 24 • p, p), where p = .001929. 

This theorem implies that the repulsion inequality fails at some density in dimension 24: 

Corollary 2. For large enough n, there exists some r so that when n spheres with centers 
xi,...x n are placed uniformly at random in T 24 , 

E [vol(U™ =1 B(xi, r)) | d(xi,Xj) > r for all i ^ j] < E {vol{lff =1 B(xi,r))] , 

where B(xi,r ) is the closed ball of radius r around the center Xi. 

In other words, conditioning on the pairwise repulsion of the centers of the spheres can 
decrease the expected volume of their union! 

Note that working in the torus is not essential to the result: the same holds if the centers 
of the spheres are chosen at random in a box in M 24 large enough so that boundary effects 
are negligible. 

Proof of Theorem [H Consider some packing of n spheres of radius r p in T d with density p. 
Place a sphere of radius (1 — t)r p for 0 < t < 1 around each center of the packing. If we 
place a new set of centers, one in each of these spheres of radius (1 — t)r p , then the spheres 
of radius tr p =: r/2 around them will be disjoint. The density of such a configuration is 
a = nvd(tr p ) d = t d p, or in other words, t = (a/p) 1 / d . We can lower bound the probability 
that n random centers of spheres of radius r/2 will be disjoint by the probability that each of 
the n centers falls into a distinct sphere of radius (1 — t)r p around the centers of the packing: 

Pr [G d (n,r) is empty] > ^((1 - t) d p) n = ^ (l - (a/p) 1/d ^j p n 

and 

— — logPr[Grf(n, r) is empty] < 1 — dlog ^1 — (a/p) 1 ^^ — logp + o(l). 

The birthday inequality, however, asserts that 

-logPr[Gd(n,r) is empty] > 2 d ~ 1 a. 

For d = 24, there is a sphere packing of M 24 of density via the Leech lattice [15, 5]. For 
any e > 0, and all large enough n, we can find a packing of T 24 with n non-overlapping 
spheres which has density at least — e - Choosing p = .001929 ~ ~ 5 • 10~‘, we can 

compare the birthday inequality lower bound on the free energy BI(p,t ) to the cell model 
upper bound CM(p,t). 

F(t ) := BI(p,t ) - CM(p,t) = ^(2 t) 2A - 1 + 241og(l - t ) + log(p). 
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A calculation gives F(. 79) > 0 and F'(t ) > 0 for t £ (.79,1). This proves Theorem [2j Corol¬ 
lary [2] follows immediately since the sequence of repulsion inequalities implies the birthday 
inequality. □ 

The hard square model is the hard sphere model under distance: configurations of 
disjoint d-dimensional axis-parallel cubes. Cubes pack particularly nicely, with a maximum 
packing density of 1. We use this to show that the birthday inequality fails in all sufficiently 
high dimensions. 

Theorem 3. The birthday inequality fails in the hard square model for some range of densities 
in dimension d > 6. For d = 6, the inequality fails for a £ (.4, .95). For d > 6, the inequality 
fails for a £ (v d ifid)’ where 77 ^ ~ 2 1 ~' / log( 2 ) • d as d —>• 00 , 77 ^ —»• 1 as d 00 . 

Proof. Here we ask for a given d if there is some a € (0,1) so that 

1 — dlog(l — a 1 ^ d ) > 2 d ~ 1 a. 

A numerical calculation for d = 6 and some calculus give the theorem. □ 

3. The hard-core model 

Assume n is such that n 1 / d is an even integer. Let Z d (ri) be the d-dimensional discrete 
torus of sidelength n 1 ^ (with a total of n sites). Assume a is such that an is an integer. We 
define a fixed-density hard-core model as follows: 

Definition 6 . The fixed-density hard-core model HCd(n,a) for a £ [0,1/2] consists of a 
uniformly chosen random independent set of k = an sites in Z d (n). 

This model is a natural discretization of the hard sphere model. It is closely related to the 
hard-core model with an activity parameter A: an independent set I C Z d (n) chosen with 
probability proportional to A^L By conditioning on |/| = an we obtain the fixed-density 
hard-core model defined above. In the terminology of statistical physics, the fixed-density 
model is the canonical ensemble, while the activity parameter model is the grand canonical 
ensemble. 

We can define the partition function and free energy of the hard-core model: 

Definition 7. The partition function, Z d (n,k), of the hard-core model is defined as: 

(6) Z d (n, k ) = IS(k) := ff of independent sets of size k in Z d (n). 

k 

We can write Z d (n,k) = Pr[A^ is an independent set] where is a (multi)-set of k 
independent and uniformly chosen sites from Z d(n). 

Definition 8 . The free energy of the hard-core model at density a is: 

(7) F d (a)= lim - log Z d (n, an). 

n —>00 n 

Note that we do not take the negative of the log partition function here, and so while we 
obtained lower bounds on the free energy of the hard sphere model, here we will obtain upper 
bounds. 
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We can write the free energy in terms of the probability that is an independent set: 
(8) Fd(a) = a — a log a + lim — log Vi[X an is an independent set] . 

n—>oo Tl 

We can also define the fixed-size hard-core model on the d-dimensional Hamming cube Qd, 
with vertex set {0, T} d and edges between vectors that differ in exactly one coordinate. The 
partition function and free energy are defined as in the hard-core model on Z^(n). 

Our first theorem of this section is that the repulsion inequality (and thus the birthday 
inequality) holds in the hard-core model at a sufficiently low density on any d-regular graph. 
We consider a rJ-regular graph G on n vertices, and select a set of k vertices Xf. uniformly 
at random with replacement. We use the convention that two vertices in X^ form an edge if 
they are neighbors in G or if they are identical; so if Xf. has no edges, it is an independent 
set of size k in G. Thus we have p = , the probability that two randomly chosen vertices 

form an edge. 

Theorem 4. For the the hard-core model on any d-regular graph G on n vertices, at densities 
a < (d + l)~ 2 , the repulsion inequality (0|) holds. 


As a corollary via the birthday inequality we get an improved upper bound on the number 
of independent sets of size an in all d-regular graphs, for a < (d + 1)~ 2 . 

Corollary 3. For a < (d + l) -2 , the number of independent sets of size an in any d-regular 
graph G satisfies: 


IS(an) < 


n 


(an)! 


1 - 


d + 1 


n 


/ ocn\ 
V 2 ) 


On the scale of the free energy, this gives: 


( 9 ) 


— log IS(an) < a 
n 


a log a — a 


2 


d+l 

2 


For a < (d + l) -2 , Corollary [3] improves the bound given by Carroll, Galvin, and Tetali in 
[3|R . Specializing to Z^(n) and Qd we get upper bounds of 

a (1 — log a — a(2d + l)/2) and a (1 — log a — a(d + l)/2) respectively on the normalized log¬ 
arithm of the number of independent sets of size an. As far as we know these are the best 
bounds known on the number of independent sets of a given size in Z d(n) and Qd at these 
densities. 


Proof of Theorem [7} The proof is essentially the same as the proof of Theorem [fl and this is 
one of the motivations of this work: to find methods for analyzing the hard-core model that 
generalize to the hard sphere model. 


Let 14 be the fraction of vertices in G at distance at most 1 to a set of k randomly chosen 
vertices. Let Ef. be the event that the set of k random vertices is at pairwise distance at 
least 2 in G. We can assume that k > 2 (and thus n > 2 (d + l) 2 ) since the case k = 1 is 
immediate. Let p = the probability that two randomly chosen vertices in G coincide or 
are neighbors. We will prove the following estimate: 


( 10 ) 


E[Vfc| E k ] > 


k(d + 1) 
n 


fk\ d{d - 1) 
\2y n 2 (l — kp ) 2 


2 For this range of a , the best upper bound in [3] on IS (an) is the third bound given in Theorem 1.6, 
2“ n ("' 2 ). On the scale of the free energy this is —alog(a) — (1/2 — a) log(l — 2a). Some calculus shows that 
the bound © is lower for a < (d+l) 2 . 
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where the randomness is in selecting the k vertices uniformly and independently at random. 
By inclusion/exclusion we have 


E[Vfc] = 1 — (1 -p) k <kp- 


p 2 + 


k\ p 3 _ + 1 ) 

3/ n 


{d + lf 


ru 


1 - 


{k — 2)(d + 1) 
3 n 


and using (fTUT) we get 


E[V k \E k ]-E[V k \ > 
> 



d(d-l) 1 

(d + l) 2 (1 — kp) 2 

d(d- 1) 1 

(d + l) 2 (1 — kp) 2 


(fc-2)(rf+l) 
3 n 

kp\ 

"3" ) ' 


This is non-negative when a = k/n < : the RHS is decreasing in k, and so it is enough 

to prove when k = n/(d + l) 2 . This follows from an elementary calculation and proves 
Theorem 01 modulo the estimate (1101) . 

To prove (flT)]) . we use inclusion/exclusion again to bound 


E[V k \E k ] >kp- Qe[F(1,2)| E k ], 


where 1/(1, 2) is the fraction of vertices in G at distance 0 or 1 of the first and second of the 

1 2 

k randomly selected vertices. Let A v ’ be the event that vertex v neighbors both of the first 
two selected vertices. We write 


E[F(1,2)| E k \ 


E[F(l,2)-l gfc ] 

Pr [E k ] 

1 Pr[A v ’ D E k ] 

71 v&G Pr [- Bfc ] 


1 ^ Pr[^i’ 2 n E 2 \ • Pr[E k \Al’ 2 n E 2 ] 


n 


veG 


Pv[E 2 }-Pi[E k \E 2 ] 


n 


J^Pr[Al’ 2 \E 2 ] 


veG 


P*[E k \Al’ 2 n E 2 ] 
Pr[E k \E 2 ] 


If the neighbors of v form a clique then that term in the sum is 0. We assume from here that 
there is at least one edge missing from the subgraph of u’s neighbors. 

Consider Pl nE 2 ] . Again we have Vx[E k \A^' 2 n E 2 \ < Pr[^_ 2 ], and 


Pi :[E k \E 2 ] > 


> 


Pr [E k -2] Pr[E k \E k - 2 \ 
Pr [E 2 ] 

Pi[E k _ 2 }(l - kp) 2 
1 ~P 


which gives 
( 11 ) 


Px[E k \Al' 2 n E 2 \ < 1 -p 

Px[E k \E 2 } - (1 -kp) 2 ' 
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Next note that for any d-regular graph G , 


( 12 ) 


i^Pr[^- 2 |E 2 ] 

veG 


1 ng)-3-#{C^inG} ^ d 
n (”) - dn/2 ~ n 


Inequalities (fill) and (fl2l) give 


K[V(l,2)\E k \ < 


d(d- 1) 
n 2 (l — kp) 2 


and thus (1TU1) . 


d- 1 
n(l — p ) 


□ 


We now show that the birthday inequality fails in general for d-regular, bipartite graphs 
with d > 6. 

Theorem 5. For d > 6, there exist constants a l d £ (0,1/2) so that for n large enough, the 
birthday inequality fails for the hard-core model on any d-regular, bipartite graph G on n 
vertices at densities a € [a l d , 1/2]. Asymptotically, a l d ~ 21og2/d as d —> 00 . 


Proof. For a lower bound on the number of independent sets of size cm in G, we use the 
parity lower bound: any subset of one side of the bipartition is an independent set, and so 

IS (an) > f™/ 2 ) . 

\ an J 

The corresponding bound on the free energy is 

Fd(a) > —a log(2a) — (1/2 — a) log(l — 2a) + o(l). 

The birthday inequality asserts the upper bound: 

Fd(a) < a (l - log a - a ^ 1 j + o(l). 


Some calculus shows that these bounds cross for d > 6, and that asymptotically as d —> 00 , 


the crossing point is a d ~ 2 log 2/d. 


□ 


4. Matchings 

In this section we use the repulsion inequality to give bounds on the number of matchings 
of size k in a d-regular graph G on n vertices. Such a graph has nd/2 edges, and each 
edge shares a vertex with 2d — 2 other edges. We let p = the probability that two 

uniformly chosen random edges (with replacement) coincide or intersect at a vertex. Then 

the birthday inequality asserts that PrfE'fc] < (1 — p)^\ where Ef, is the event that k edges 
chosen uniformly at random from G form a matching of size k. The repulsion inequality 
states that E[I4|£/] > E[V%], where 14 is the fraction of edges covered by a set of k edges: 
the fraction of edges that are contained in or intersect the set. Since a matching in G is an 
independent set in the line graph of G, we could apply Theorem [3] to get a bound, but we 
can do better working directly, since for d > 3, the line graph of a d-regular graph contains 
many triangles. Let M{k) be the number of matchings consisting of k edges in G. We show: 
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Theorem 6. For a < the repulsion inequality holds for matchings of size aff in a d- 
regidar graph on n vertices, and as a consequence 


M(an/ 2) < 


(nd/2) an l 2 
(an/ 2)! 



(cx.nl 2 \ 

2d -i y ^ ) 

nd/2 ) 


On a logarithmic scale, this gives 


(13) 


— log M(an/2) < alogd — a log a + a 


a 2 2d — l 
~2 d~ 


For a = 0(d 1 / 3 ), Theorem [U] improves the bound given by Ilinca and Kahn in Efl 
Together Theorem [6] and [11] show that the birthday inequality holds for matchings of all 
sizes in d-regular graphs. 


Corollary 4. The birthday inequality holds for matchings of size k, for all k, in every d- 
regular graph on n vertices. 


In would be nice to prove that in fact the repulsion inequality holds for matchings of any 
size in a d-regular graph. 


Proof of Theorem 0 Let m = nd/2 be the number of edges of G, and p = 2d m l . We want to 
show that E[Vfc|Kfc] > E[14]. Let L 2 be the number of edges of G that are covered by two 
edges of a matching of size k but are not part of the matching themselves. Then 

E[14| E k \ = k + (2d - 2)k - E[L 2 \E k \ = m (pk - E[L 2 /m\E k }) . 


By inclusion/exclusion we have 


E[Vfc] = m(l — (1 — p) ) < mpk — mp ( ) + mp 6 


< m [ pk — )p z [ l — 


kp 


So it is enough to show that 


We can write 


kp 


E[L 2 |K fc ] < m p 1 - — 


E[L 2 \E k ] = J2 ( 1) Pr[4’ 2 \Ek] 

e£G ^ ' 


3 The bound in m, translated to natural logarithms, is 
4 log M(an/2) < a log d — a log a — 2(1 — a) log(l — a) — a + (logd) /(d— 1). Subtracting the first two matching 
terms then power expanding around a = 0 gives a — a 2 + (log d)/(d — 1) — a 3 /3 — ... for the Ilinca-Kahn 
bound and a — a 2 + a 2 /(2d) for the birthday inequality bound (1131) . These cross when a = 0(((log d)/d) 1 ^ 3 ). 
In particular, for all larger a, the Ilinca-Kahn bound is stronger than the birthday inequality, giving Corollary 

m 
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1 2 

where A e ’ is the event that edge e is covered by edges 1 and 2. Now it is enough to show 
that Pi/A^C j-E/] < p 2 ( 1 — kp/ 3). As in the proofs above we write 


Pi[A 1 /\E k ]=Pv[A 1 e ’ 2 \E 2 \ 

< Pt[AI’ 2 \E 2 ] 


Pi-[E k \Al ’ 2 n e 2 ] 
Pr[E k \E 2 ] 

1 — p 

(1 — kp) 2 ' 


We calculate 


Pv[A\’ 2 \E 2 ] 


2{d — l) 2 

m 2 { 1 — p) 


which gives 


Pr[A l f 2 \E k ] < 


2{d-lf 


m 2 ( 1 — kp) 2 

Our assumption that a < 3/28 implies that kp < 


A, and so 


Pr ^' 2 ' E ‘i £ 


2 — 


(2d - l) 2 


777T 


1 _kp 


which shows that the repulsion inequality holds. 


= p 2 (l — kp/ 3) 


□ 


5. Conclusions and conjectures 

We conjecture that the lower bounds on the density at which the birthday inequality holds 
in Theorem [T] and Theorem [I] can be extended to the entire fluid phase of the hard sphere 
and hard-core models. 

We describe two notions of the fluid phase of the hard sphere and fixed-size hard-core 
model. The first is decay of correlations: 

Definition 9. Let xq, x' 0 , xt be positions in T d or lattice sites on Z^(n). Let Ao (resp. 
A' 0 ,A t ) be the event that the position xq (x' 0 ,xt) is covered by a sphere in the hard sphere 
model or occupied by a particle in the hard-core model. Then the model exhibits decay of 
correlations at density a if there is some constant c a > 0 so that 

|Pr[Aj|A 0 ] - Pr[A f |A q] | < g(c a d t /r) 

where dt = mm{d(xo,xt),d(x' 0 ,xt)} and g(s) is some function so that linn^oo g(s) = 0. 
(For the hard-core model, we take r = 1). The model exhibits exponentially fast decay of 
correlations if we can take g(s) < e~ cs for some c > 0. 

The second notion is the rapid mixing of a specific Markov chain with the hard sphere or 
hard-core distribution as its stationary distribution. One such chain is the single-particle, 
global-move dynamics (see e.g. M). A single move of the Markov chain consists of selecting 
one center of a sphere or one particle on the lattice uniformly at random, then selecting a 
position or a site uniformly at random from T d or Z ^(n) and moving the center or the particle 
to the new position as long as it does not violate the hard constraints of the model. We say 
the chain mixes rapidly if the mixing time is a polynomial in n. 

Conjecture 1. If the hard sphere or hard-core model is in the fluid phase (say it exhibits 
exponentially fast decay of correlations or the Markov chain above mixes rapidly), then the 
repulsion inequality (and thus the birthday inequality) holds. 
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The intuition behind Conjecture [T] is that at a sufficiently low density, conditioning on 
particles being repulsed from each other should have an essentially local effect, and locally, 
conditioning on repulsion increases the volume covered by the union of balls around the 
particles. However, beyond the fluid/solid phase transition, long-range correlations come 
into play, and conditioning on the repulsion of particles can force them into global lattice-like 
configurations with holes, and thus the volume covered may actually decrease. Note that a 
model of random matchings of a given size on the d-dimensional lattice, the monomer-dimer 
model, does not exhibit a phase transition [TO] , and Corollary [I] shows that the birthday 
inequality holds at all densities. 

Conjecture [1] has several consequences. First, it would give a mathematical proof that 
the hard sphere model in dimension 24 undergoes a phase transition, which to the best 
of our knowledge has not been proved yet in any dimension. The density at which the 
birthday inequality fails would mark an upper bound on the critical density for the model, 
as exponential decay of correlations (or fast mixing) could not hold. 

Second, it would imply that the critical density in the fixed-size hard-core model is upper 
bounded by a l d from Theorem 0 in particular showing that a phase transition occurs at 
densities 0{l/d) on Z d (n). The best known analogous bounds in the hard-core model with 
fugacity parameter A are A c = 0(d _1//3 ) given by Peled and Samotij m improving the bound 
of 0(d _1,/4 ) from Galvin and Kahn [8]. Proving Conjectured] would give the optimal bound 
up to a constant factor A c = 0(d _1 ) as the typical particle density a in the hard-core model 
with fugacity A is bounded below by 

The next conjecture concerns independent sets in d-regular graphs. 

Conjecture 2. Suppose 2d divides n. Let H dn be the graph consisting of n/{2d) disjoint 
copies of K d d , the complete bipartite graph on two sets of d vertices. Then for all 1 < k < n, 
H d , n minimizes the expected number of neighbors of a uniformly random independent set of 
size k over all d-regidar graphs on n vertices. 

In the notation above, Conjecture [2] states E# d n [I4|£T] < for any d-regular G on 

n vertices. Conjecture [2] immediately implies a theorem of Kahn [12] and Zhao [22] that H d n 
maximizes the total number of independent sets in any d-regular graph on n vertices, and 
implies the conjecture of Kahn [12] that H d n in fact maximizes the number of independent 
sets of size k, for every k. Conjecture [2] is in fact much stronger: it states that H d)U maximizes 
the ratio of the number of independent sets of size k to the number of independent sets of 
size k — 1 for all k. The elementary proof of Corollary 0 suggests that proving Conjectured] 
may in fact be easier than trying to bound the number of independent sets of a given size 
directly. 
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